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A53Ta.lCT 


Resistance to flow on a ri?id flat beds having 
randomly spaced unifoim sand grains and densely packed 
nonunifoCT sand grains was investigated, to cover the 
transitional state of flow. The median size of the 
sand grains used was of 1.5 mm. The roughness of sand 
beds are represented in terms of Nikuradse’s equivalent 
sand grain roughness Kg/d^^ and ^ift in velocity scale 
ji.u - In d^o • Based on the present experimental 

results and results of previous investigations carried 
out at this Institute, models to represent the rou^mess 
of sand beds in terms of geometry of the beds represented 
by rou^njss concentration and nonunifonalty parameter 
«r/d^ and flow parameter represented by Grain Shear 

ReTOolds npnbor — r has been developed. These 

v 

models predict the roughness of sand beds f airely 
well. 
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3j^, Gonstnoits 

Man tar of the sand ^ain 

Brain sizse through #iich 50 percent of saople is passed 

Grain size throu^ i^ich 65 percent of aaaple is passed 

Grain size through which 85 percent of sanple is passed 

Grain size through vtoich 90 percent of sanple is passed 

Meaa height of roughness grains 

Jikuredse' s equivalent sand grain roughness factor 

rime average velocity of air at any height y 

i?ree stream velocity of air in the test section 

,To , 1/2 

Shear velocity, 7* = 


Shift in velocity scale 

Grain Shear Re’/nolds nun her based on ^50^ Be* = 

P ?rpendicular height from the theoretical bed level 
(y^ = y-G) 

Distance (Ordinate) measured perpendicular to the 
bed from the top o*" the grain. 

Zeroth man ant of spectral function 

Second moment of spac'ral function 

fourth moment of i^ectral function 

Spectral parameter 

Significant height (Hg = 4- j^) 



a v/ara-je 3tnplitude of sand ^rain beds 

(j~^ ‘Standard A-rrlation of amplitudes 

s A n rsi^Q sp ac i n« 

.'itraodard deviation of spaein^s 
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y: 

Pa 

CJ 

Yf 

Ys 


3oundarv layer thickness mesisured trom. the top of ^ain 

Boundary layer thickness measured from the theoretical, 
bed level (6.^.- 6-6) 

The location of the apparent origin for velocity 
distribution abo-j'e the smooth surface of the flat bed 

Von Karman's universal constant 

Mass density of air 

Standard deviation for the dinmeter (size)of sand 
grains 

Kinematic viscosity (air) 

Roughness concentration 

Specific weight of fluid (air) 

Specific weight of solid (sand grains) 
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A, Introduction 

?he flow phenomena in rivers, canals and deserts 
are controlled by roughness of sand beds. Before initiation 
of motion of sand grains, roughness of sand beds can be 
considered similar to flow over rigid beds, dikuradse, as 
early as 1953, developed methodology for rough wall flow 
analysis by investigating flow in pipes having their walls 
roughned with unifoiTT sand grains glued closely. Schlichting 
developed roughness scale for other roughness which are 
less densely placed and called it as equivalent sand grain 
roughness and is denoted as 5g, Th^re was no systanatic 
investigation on the roughness of sand beds till recently, 
Investigat ons carried out byP.K, Xittal and 3. Mittal 
showed the importance of nonimi-^ornity of sand trains 
in sand bads, for completely rough turbulent flows. David 
investigated the rou^-ness of sand beds for transitional 
and sraooxh turbulent flows. It was found that roughness, 
denoted in terms of Nikuradse's equivalent sand grain 
roughness varies with state of flow and aonunif ormity 
in sand grain size distribution, fhere was wide gap in ths 
transition region to fully rough ttnrbulent regions in the 
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aoove exp>3rim jnt:al series, H--;re an artonpi to investigate 
this reT'on has bee a made. 

•^''Ughness of san'l beds can be predicted by using 
CO 'relations curves davjlopijd by Da'/id, An attempt to 
develop model for rough turbulent ■''lovu on sand bed vras 
made by David. This model can not be extended for 
transitional state of flow. Development of model to 
predict roughness of sand beds for rough turbulent aM 
transitional regions was of necessitv, H 'ro as attempt 
to develop model for rouj^ness scales in terms of geometry 
of bed surface and flow properties has b ^on att;aarted. 

B. jlelavant Liter^.t!ure R .view 

Jikuradse developed a fr-ane work of rough wall 
flow analysis by investigating flow in sand roughned pipos. 
liikuradse’s investigations showed that the affect of 
roughness on the velocity profiL* was confined to a thin 
layer near the wall. *^0 velocitv distribution in this 
layjr follows a logarithmic •‘"unction. This is universal 
in character and is determi’aed only from -the wall conditions 
and distance from the wall, According to the current ■view, 
the velocity distribution in the logarittmic zone does not 
differ \diathcr it is boundary layor flow, pipe flow or free 
surface flow (Monin and Yaglom, 1971; Hinze, 1975). 
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Jr i.:ural3'j’3 invisti rations, th-i ''r.c of 
Schlicht?.n--:'s (1936) exp ^rim j.its rjvoal^d that , for 
rourh ourfacea h-' 'int Ijssjr o jncc.‘?i'tration el. 
the rou,dio.jS3 parrsa ■•‘•ir was '’■‘''’ereit ■f'ron th ,• jjir * o*? 
th': roag*'.ajs?. jljr jit. Jc Jlicat 3f o -td that thj 
roughness effect is a function of th5 average height of 
the rou^ness elaa nts and also of their site, shape and 
concentration distribution over the wall surface, therefore, 
in order to characterise rouf^ness, hj made use of 
' 'Jikurads j equivalent sand grain roughnjss, Kg’* viiiich is 
defined as the size of the uniform sand grjdn that, produces 
the same resistance coefficient as the actual roughness under 
the same •^’low condition. Using Kg, to ropr :sent both 
unif 02311 and nonunifoim irougnness, the velocity distribution 
in thi- logarithmic zone is written as ‘ 

(yAs) + B (f, Kg/si ) <i) 

Based on the concept of Nikuradse, to rtpresont Kg in 
toras of particular 3ize of send grains, is usual practice 
in sadim nt transport stuiies, To oniao a few of the studios; 

Hinst^in and Barbaixjssa (1952) considered Kg = dgg 
Simons onfl Richardson (l%6) utilized KgsdLgg idicrcas 
Kamphuis (1974) used Kg « Sd^Q, 
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'Uil brook ’"id (l’"*?) lirr gtiT^t-jr low 

rysi-taicj ’.wall be affected V-^' addin'? a stall bu*^ d 'finite 
p ’’opo rti 3’o. od I'T^e grains to unifor’^ ■fine tr-Jl?is. 'ffhov 
ob^er --id db-'d the p-''?3Mce O'" l'ir'2r Tr‘'’i.n3, "bouTh only 
in aiic-ll proportion, consit jr nly off--fC*£i the roriot' ic 3 
mech .niat cue to th'ir shielJi ir ed ’oct on fin re oroiao 
and thus reduces the t!" jctlvoness of atelier ’reins. 

Iheir res'ults ’tlso indicate that rou?hn;ss conceit ration 
is not 1 ?ood para'3?t3r to deal -.vith nonuni'fom rourhaoss, 
7/hile it is us3.ful for unifona rou?ha-;es. IPhis is b'cause 
'cone-, itr :t ion' bocornss vatac Then the whole rsa o,: bed 
is covered bv s.nd rraina. 

O'Lou'.ilin anJ 1-Iac tonild (ib64) did experim .ts 
on said ^p-iins os routhn'osp '1 'to "its randomly spaced 
in th ; var.ljus o onoantrations. .^ccordin? to tti.vfl, the 
di^-’eronce in resi=’'‘'ince □""crad by ro^ulirly arranted 
spher-a and irre'Tulorlv arr-iniecl san.i grains is due to 
the irrb palarit:?" in sh'-'p.e and the raniomness in pTr^ngenent. 

It mav bo noted that ropres ‘ntation 6f Kg in 
terms of xiy pa-rticular size is not proper, itisteatl 
pframi^t^rs ieseribin? the c istribution of sand prains 
namely median size d^Q and standard diviation^Sr^^'^'’ ^3 
better par«*i 3 t rs.* P.K, Mittal ^pd S;. Ilitrtal conducted 
experiai'ents to find th'3 effect of nonuni^oimity in 
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■?r?in siTie repr'iss-itei in terns of c5'/a^Q on. the r^urhnesa 

scale Kg. The 7 observed tliat for coarse sazid grains , ^g/^'j0 

is a onlaue ■^'onction o' nonuni^'omitv coefficient 

and is iniepenient of Ct-nin Shear "lo -nolds i~her, - 

David conducted a throu.^ inveeti nation to find -^he 

den '♦ 

effect of Srain Shear Revnolds nunber — ^-rr — alon,? 

""p 

^^50 ^^s/^50' ^avid sho’wed from his experinent 



for 7* den/^ = 60, 24 and 8 that the functional form 

V* dcQ 

Kg/d^o with ®/d 5 Q donot coincide with higher —n^—— 
value, but follow different curves for each range. 


It may be observed that the gap in d,^ values 
between l.K. ilittal and Sudhir /Ii"tal data and David's 


data was wide, so a necessitv was felt to investigate 
for d^Q values that lies bet’vee.i 200 and 60. Also there 
is a necessity to develop an. analytical relationship 
between Xg/d^Q and ^d^Q for various d^Q . Here an 
attempt is made to investigate experim.e rtally the ef'ect 
of nonuni.formity in grain siae on rou.ghnes3 scale for 
di^o order of 100, Based on experimental data, a method 
to predict ^[g/d^Q by knowing (CT/d^o and d^Q is also 
att<anpted. 


Investigations of Schlichting (1936,)» O'Loughlin 
and .lacDoncOd (1979), David (1980) and Sarin (1980) 
showjed that the arrangement pattern, she|>e and relative 
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size of the Trough.iess element are inpor^Tiit. fhe 
roughness concentrat on defined as the ratio of the 
projected area of the grains *0 the tot j 1 area of 
the bed is used to repnesiit on* ot* the geometrical 
parameters of iLiiforra rjughness eln-’iits ft.en sp-irc:*!/ 
distributed on the bed. In the case of densely packed 
sand grain beds, the nonunif ormity paran^t ?r *5" /d^Q 
is considered to represent the geometry scale of the 
sand bed , 

G. Shift in Velocity and Length Scales 

Another form O'*' representing roughness, 5.3 in 
terms of shift in velocity scale and shift in length 
scale. 


This method was intro iuced by Kama (1954). Here, 
the roughness ef-*’ect is considered equivalent to shift 
in the velocity profile from smooth wall to rou^ wall 

by a value du/Y* = and is written as 

u . n A u 

( ) - ( ) _ ( 2 ) 

'/♦rough '■ y '^smooth ' ' 


where 


and 


u 

(r ) smooth 
^ # 


u 


('^-)rough 
*+ J 

in which y = ~rT" 


: ^ In y+ + 3^ 

~ In yAs % 

• ; 3g and \ are constants for smooth 
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anc conpl=;t^ly rough beds. 7sln? th?sa 3 tp”*?S 3 i-’' is, 
the shift in vsluoitv may be represe.ited in terms of 
33 and as 

Au"^ = In Xg + (Ig - 3t^) ( 4 ) 

Hama related Kg with u^ as 

£xn^ = 5.6 log (K^ + 3 . 3 ) - 2.92 

o 

?or fully rough turbulent flow condition, Hama 
( 1954 ) and ^lauser (1956) detenined the value of u'*' 
for different types of roughness, latterman showed that 
u'*’ - ~ In is a f unction of roughness concentration 
for two-dimensional rohghaesses. 

Perry et al, (1969) investigated the following 
relation 

1 V* e 
Au+ = — 

The value of 3 was foun5 to be varying with 
geometry of the roughness scali and =»l30 with state of 
flow (d 5 Q ). 

P.K. iittal, S. Hitt el and David showed that the 
function In K"*" is a function of and 

for nonuniform sand grains densely pacified, David, 
Sarin and Aslaoi showed that the function u — In ^50 
Is a function of rou-^bness concentration and d 5 o for 
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unifora trains ri.xiTal;^ spaca-i (actio? as 3-iiKi 30 s iooal 
ro ughnessas) . 

D. Theoretical lad Lavel 

The y-looity distribution on a roufh wall bchavas 
as if its ori'?in is loca+ed at some distance, 3^ below 
the crest of the roughness aleraents ( Moore, 1951), 

Blineo and Parthenir'es (1971) found a relation for 
uniform densely packed bed 

3^ = 0,21 Sg (7) 

Einstein and El-oamni (1949) had also obtained this 
relation for pebbles. Kaaphuis (1974) assumed that the 
virtual bottom is 0.7 d^Q above the plane to which the 
sand grains were attached. David showed that the S/d^Q 
is a function of ^^/d^Q and d^Q (\i4icre '9' is a distance 
measured from plane surface to the theoretical bod level). 
For d^Q values greater thaii 220, the value of 9/dgQ was 
less than unity. Tor d^Q < 220, it was observed that the 
theoretical oed level was much above d^Q level. David, 
Sarin and Aslan showed that the theoretical bed level 
coincide* with the geometric bed level for uniform sand 
grains and glass beads sparcelv distributed for rou^ 
turbulent flow conditions. 
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E. Present InveatigatiOR. 

An investigation on the roughness chiracteristios 
of sand grain beds was planned as a continuing study. 

The nonuniformity parameter defined as 5” 7^50 varied 
to study the rougjmess characteristics with nonunifoiBiity 
of sand grains. 

Under this prograome, P.K. Mittal (1977) carried 
out experiments on rough surfaces represented by sand 
grain beds having a median diameter (d^Q) of 1 mm, 2 mm, 

4 mm, 6 mm 8 mm. The nonuniformity of bed was obtained 
by varying the standard deviation of sizes to 2-3 values, 

S. Mittal (1978) worked m One median size 

cr- 

SI 5,0 mm ) but varied to about 10 different values. 
David (1980) investigated 3 series of sand beds having 
d^o = C*14 mm, 0,39 mm and 0.925 mm. Each series had a 
number of beds with different Sarin investigated 

the effect of rou^ness concentration of size 0.925 mia 
sand grains and 3.0 mm glatss beads on the rougjiness scales 

Ks/d5o, -du’^ “ ~ *^50 ^ ®/^50* 

The present work is planned to investigate the 
following aspects: 

(l) The effect of laanunif oisnity in sand grains 

size with median value 1.50 mm to relate the 
work of P.K. Mittal, S. Mittal and Bavid. 



Bxperiments of uniform size vsdLth d^Q = 1,5 nm 
for different roughness concentrations v/ere 
planned to find the effect of randomness ia 
spacing on roughness scales. 

To develop a geometrical peuram jter to represent 
the roughness scales for above oases. 

To develop a method to predict Kg/di^Q by 
knowing and d^Q , 
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A. Details of '^ind Tunnel 

An open circuit wind tunnel vrith closed test 
tube' section is used, A honeycomb is located at one end 
of the wind tumel along with five screens of size 
1,25 m X 1.25 m throu^ #iich air enters into the test 
section, The hoheycc3inb-cum-screen portion is 2.3 m long, 
after wdiich there is a contraction which leads to the 
test section. The test section is 0,4 m square and 4 la 
long. Prom the test section there is a gradual expansion 
which leads to a circular section of diameter 1,0 m, At 
this end is fitted an exhaust fan of 1,0 m dia wtoich 
sucks air throu^ the test section. The exhaust fan is 
operated by a motor of 3.7 kw, 5 hp and 960 xpca. 

For controlling the velocity of air in the test 
section, there is a pair of adjacent gates fitted at the 
end of the test section, 'fhen the gates are completely 
closed the maximum possible velocity is obtained. The 
gates can be operated and kept at a desired opening to 
get four different velociHea aispnoxiaiatelr'IS.O m/seo, , 
17.0 m/seo,, 14.0 m/sec and 11,0 m/sec. 
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J'iv ‘2 stations in thi isst section it nistauicas 
of U,14- m, 1.04 n, 1,94 ti, 2,84 n .ml 3 •73 n firom the 
le ii'in? ei o-*’ the test section were chosen at v^ich 
rteasavenj-it of vslocitv was t ilvan. s-.dll hole Is 
proYilei on the sife of the tost section at ^ac.^ of 
ttese sections for static pressure noasuren'-nta* A 
total head tube, naving. 9 nin outside diameter and 5.0 OB 
Ion? is used :®'or ths measurement of total head pressure. 
This tote.l head tube is attached to a traversin’ 
m ^clianian w’nich can move in the vertical cJ.rection with 
0.05 mt least count. The total head tube and the 
static pressure tap hole on the side the test section 
are connected to a differential mmora?-*' *r which can read, 
to an accuracy of 0.00254 cm o" wat 3r. .i layout of the 
wind tunnel is shown in ?i 1. 

3 • Preparation of the Rouda Bads 

i set of two plane and smooth wooden boards, 
compl-'tely occupving the test section floor of the wind 
tunnel, were used to prepare one rou.ih bed. A thin coat 
of p ant was evenly applied over one side the board. 
Roughness elpmrnts of known vs^ight. were sprinkled over 
this surface to ^t a random arrarr=tement« The bed was 
allowed to dry, without disturbance, for about 40 hrs. 
When the paint dried coraple'^aly, the bed was upturned 
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se ~he ou^. Ch*? \/eii?ht 

o* th-5 e’^’.o’cin? vere iis rl button 

cur'/Qs are riven in. i'ig, 2 ani Pla^e 1. CJains ■‘'h.ese 
curves, the vslue of nonuni^ sriity* oae.ffici-?nt <^^50 
was calculate'^-, i’or beds having unifonr grains r^-idonily 
spaced, the spacing was aeasured. .The distribution of 
sp.acing values for each bed is slWTO in I^ig. 5 and 
Plate 2. 

C , Meas ureci^nit of /elocity Praflles 

fhe 'wsoden boards constlt uting the rough beds 
v/ere laid on the test section of the wind tunnel and were 
fixed to the floor of the te«t section by 0 set of bolts 
which were flush with -^he painted surface of the bed, 

The leading edge of the bed a th taper in or3er 
to 3 oin r/ith the floor o:' the tunnel. 

'^he desire'^ free s^rea-^ velocity throng the 
test section was obtained by adjusting the gates to a 
particular ope.iing. t'our differ mt free strev^r velocities 
,were used for each bed, "which were app'^xiJiaT jly 19. 0 a/sec., 
17.0 m/s3c,, 14.0 la/sec iid 11,0 m/sec, ?or each of these 

free stream velocities a minlaua of one veloeity profile 
was taken along the centre line of the bed at aach of the 
five stations. On the average, the nucabor of velocity 
pro'^iles Tjakan on each bed ranged between. 16 to 20, 









19 


^or velocity noaauTera ^nt , t’jrj lanora .t ’ivs (dapendiiig 
upon the ranja of presauro) end a total h?ad tube (or probe) 
vd-th a verniar scale and slow r»ovin? divice was used. The 
static pressure end of tbo nanometer was connected to the 
static tap hole on the si'^n of the test section and the 
other end was oomected to the total h'ad tube. The room 
temperature was noted for each test. 

Once the req.uired free stream velocity was mainx idned 
in the test section, the probe was lowered down to the bed at 
each station. ];Iear to the bed, the pressure measurements 
were taken at very close intervaQ.s (0,5 ton or less), and 
at increasing distances from the bed. The interval was 
gradually increased to 1 mm, 2 mm, 3 nio, 5 md etc. since 
that part of the boundary layer covering the wall region is 
very important in data reduction, c ^ro was taken to get a 
minimian of 10 readings in the region extending from 5 percent 
to 15 percent of the boundary layer thickness. 

D. M easurgnent of the Grain Surface Pr ofiles of the Sand Beds 

Two specimens oi each 'md of 1.5.B sand bed series 
had been prepared on glass plates for grain surface 
measureiiK'iits. The glass plates were of size 3 ora x 20 era 
and were given a thin coat of pain* and laid side by side 
with the painted wooden boaids used in the preparation of the 
bedsi The sand was sprinkled over the glass plates and the 
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boaaiSB at the same time and in the aaae manner that 
there was no difference between the glass plate portion and 
wooden board portion of the bed except that the wooden 
board portion was tested in thi wind tunnel while the 
glass plate portion was used for grain surface measurements. 

A sand auitface profile meter developcid by David 

was used in the profile measurements* It ossentiaXly 
consisted of a needle point which could move both in 
horizontal and vertical directions by moans of horizontal 
and vertical transverse mechanism. In both directions 
verniers could be read, the least count of the vertical 
vernier was 0.0$ mm and that of the horizontal one was 
0.1 mm, She glass plate model of the bed was placed under- 
neath the needle point. The horizontal alignment of the 
glass plates ajad the base of the apparatus could be chocked 
by a spirit level. 

The Initial readings were token by lowering the 
needle to touch the msooth portions of the gloss plates* 

Now the needle was raised ijp, traversed along the centre 
line of the bed, lowered to touch the grain airfoce at 
every 2 Bm horizontal interval apod the vertical vemler 
was read. The bottom level of the grains being known, 
the height of the grains oould bo calculatod* 
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boaarfis at the same time and in the same manner that 
there was no difference between the §lass plats portion and 
wooden board portion of the bed except that the wooden 
board portion was tested in the wind tutmol while the 
glass plate portion wais used for grain surface measurements. 

A sand staff ace profile meter developfjd by David 

was used in the profile maasuronoiits. It essentially 
consisted of a needle point which could move both in 
horizontal and vertical directions by moans of horizontal 
and vertical transverse mechanism. In both directions 
verniers could he read, the least count of the vortical 
vernier was O.OJ mm and that of the horizontal one was 
0.1 mm. She glass plate model of the bed was placed under- 
neath the needle point. The horizontal aligrxaent of the 
glass plates ajid the base of the .app'iratus could be chocked 
by a spirit level. 

The Initial recdii^s were taken by lowering the 
needle to touch the smooth portions of the ^ass plates, 

Now the needle was raised up, traversed along the centre 
line of the bed, lowered to touch the grain surface at 
every 2 nan horizontal intervsil and the vertical vernier 
was read. She bottom level of the grains being knom, 
the height of the geaina could be calculated. 
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B . Go mpil a tion O'? Datd 
a. i<Iean veloclTiy data 

‘Velocity was raeasar^d on rou-rh b'ids alonr th-? 
centre line of the test section at five statiois, for 
four velocities having a magnitude equal to 11, 14, 17 
and m/sec. approximately, -he d’mamic head in t:hfi 
differential manometer was measured as head of water 
in inches. It lias been converted into mean velocity u , 
in m/sec. Ihe velocity profiles ware used in deduction 
of wall low parameters i.e. Jhu/T^^ and S. 

b . Ded uction of wall law paran»>ter3. 

Boundary layer has b.jen devided into two regions: 

( 1 ) 'tall law region, and 

(2) Defect law region 

The wall law region has been defincjd between 
to 0.2 and rest region has been defined as defect law region. 
She wall law region from Yt/^t = 0-2 is represented 

by logsG-itfiiBic T3loci1*y 

Acoordli^ to E®a (195* ), the wall l,w relationship 

in logarithmic region has bean written as 

= _L In - t 4^-9 - dt/f. (8) 
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In the above eq'iation u and y are measured 
quantities and the value ofx, is talcen as 0,41» according 
to Haraa. Ihe shear velocity 7#, shift in velocity scale 
/;u and shift in length scale C are unknowi quantities, 
lo evaluate these three unlcuo’fln qu-'intlttes, least square 
method assisted by Ifewton-laphson procedure was used. 

From computed values of au/7* f the values of Kg /^50 



The computed data is given in Table 1 and 
other details of mean of heists, standard deviation of 
heights and spectrxaa momenta are given in ippendix. 
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lABlZl 1 


*5 . ifo , 


* Rj;s%TS 


^50 

0.05 

0.10 

0.20 

0.35 

0.65 

1 i 25 

1.75 

4.80 


S.Eo, 


0.05 

0.15 

0.26 

0.40 

0.70 




77.85 

103.96 

112.05 

133.52 

133.62 

86.10 

77.82 

72.00 




110.84 

110.26 

ivr4.-i7 

107.09 

90.35 


^50 

Au - J-li dcA 
X - 50 

e 

^50 

< T"Ka 

*^50 

1*65 

- 1.5 

1.730 

0.159 

2*45 

- 0.5 

1.400 

0. 170 

3.40 

0.3 

0.960 

0.090 

4.20 

0,7 

0.915 

0.120 

4.65 

1-0 

0.905 

0.200 

4.25 

0.7 

0.900 

0.080 

3.70 

0 i 2 

0.750 

0.110 

2.05 

- 1.3 

0.452 

0.140 


- i-in dc^ 
d5o 50 


1.05 

2.00 

2.40 

2.33 

1.72 


- 2.90 

- 1.02 

- 0.67 

— 0,74 

- 1.40 


0.690 

0.850 

1.060 

1.010 

1.730 


^50 

0.060 

0.100 

0.090 

0.150 

0.120 


5 



CiE^l'SR III 


r-^^YSIS m IT^YLCPY'lTi.' 0^ !.;C15LS 

A. General 

His chs.ptsr is to ■’.ire results 

deduoad for .turMlent flows over all rouglJJ'ss beds tested 
in tbo present investig'-tion. whenever req.'uir‘=‘d, the 
results stated by Schlichting (1?36), 0*lioughlin et al. (1964)? 
David (1380), P.K. Mittal (1977), S. Mittal (1978), S-irin 
(1980) end ^Islam (1981) ere reproduced. tfeing these data, roughness 
parameters have besfi. studled.Xofieda for roughn^es paraoetors 
(Ks/%0^ .and >- 4 ^^^0- or X '®d 

(i Q h)?.ve be?!i developed. 

B. !ie an Velocity Distributions 

Mean velocity profil'^s in a boundnzy layer flow can 

b 3 r6 p?’ e 3 n s v d as 



The velocity distribution near tha wall is 

>*• 
u 


= f( T, y/V ) 


( 11 ) 
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Th? vlociGj di : '.o-’ib ’•t’.o--: in t'-s ^ons is 


= f( y/6) (12) 

V* 

fuT'.c tio'OdL r--D '.’ ’£i"it:r uion? been doocrl'bsd 

as follows: 


(i) Prandtl^s w?JLl 1';Y; 





according to "mna 


(3 = 4- -3 for snoo -'b bi d 

( 1 ?) 


( ii) v-locisy dof^ot law 

= .Vln y/6 + B* (H) 

1'* ■* 

(g = 2,3 .Cor 3 D 0 C uh bed according to ninns), 

YTh 3 rc.i»,s B aiid 3* ar' functions of flow and 
goorcstry of rouy'nji^ss# 

':T-.>,a,p valrcity piW'^'lles obcalnod fr*')ni ‘.11 roughnsss 
beds under present tnv! selg'tlon can be represented 

in two -ays. V- loclty profiles h.-vt b-en presented only 
for station 5 e-nd for velocity of lEagnltudfc 19*0 n/® oRlT- 
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'• ’’ of. 

/ • \ u V* ji, 

\ i; o«? function of — —Ji 

Toloci p.Tofll''®. !■> t‘^r^.3 of H«. haT?" b' 'n 
V* y+ 

plott-^fi -p-,- wh<*3» 3R| “ j^3 for all the bads 

in 5’tgs. 4 - ad 5« Here Talooitj 4^®^“lb«tlo!i for “ach 
h-’d follo’?s i gt: ^ight line parallel to oa^ ■ootbex* 

Th'i v''lo2lty profilf for roai^''08 b«f in Sh'' rsgion 
of >s =.25 or ~ .65 h-»8 ths highest *'valu« of 

5'or vmi'.tioa of ^ and '*lth«^r side of theg® . 

r ■lu'^s, r=;30'.e-tl7^.1«‘, d*scr?;.u?»g ays tic 3 Ilf. The 


Tolocit^^ proftl' 3 iidiC'-'te tb.t 


l3 fUactiOTi of 


s ’Z-** ' o ^ c X 


las ^ or °€ 


( ii) as -A functi''Ti of 

Y* ^ 

..‘h' plo' a of lot hare hr-»r, 

plotted in 51g. 6 and 7* fh^se prof 11 s f vll on .» single 

tarra. Th??:a is a soatter ne'X tfc* well, vhioh is tbe 

effect of flgeoeitj and notwinifomiitf In gr-itn si 88. 

Ihs roa#d3n«as scales ^ aad S^d^are f bn o t j j aMi of 

%0 

rough.!.- AS g«oa*trf , lll?c ». gr^'ln she^r 

Refolds nai?ib«^r .r:!..-.— * fhe analfsis of th^'s# fttnotloa^l 

V 




FIG. 4 THE LAW OF WALL FOR 

















fo::?TTs is 


up in f'i 






( 1 ) 


a fu.iotlon 7 fc /6 

-rj ' t' “ I 


^h" ^=*1001^7 def'^ct 


U >“ u 

oo 


'h 


again 7 


plotted in ?iga.8 nod 9 . The data fellow's 


.7^/6 j. h.ig b’en 
reasoaablv .. 


single curve for y^/6^ > 0.2. This cl-^ra.iy indicates 
that the effset of w-u.1 roughness is obseit In th'' defCNSt 


Ibw region. In th' regi-^n y./Sj. < 0.2, the curves follow 
different curves h wing the mag^ituie of ti* dlffer'-nt 
for e soh case. This is con'>lc<ercd due to effect of 
geometry of t^- > rough b-^d eurf->c?. 


5 ^ Sand . 

•!S\1 

The roughness sc.'les K- ?jnd~— that were Intro - 

^ V* 


ducod in the velocity profiles were f'^tird ^:o vary 
with th;-' surfeoo chrtacteri''tics X or 'f/d^Q 
"lid flow chera:teristic3 


V 


Roughness p ra^eters fo^ randomly arranged, 
unifo^. or nonuniform siz' roughness clEiaenttii can be 
r^l-jsd fuactionally to flow >roperti •■-. and rough surf tcc 


prop nr tins as 


IJ 1 ^AMPUt< 
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FIG 8 VELOCITY DEFECT LAW FOR SAND GRAIN BEDS 
(V5 DS) 



5^ «i 



FIG- 9 VEIjOCITY DEFECT LAW FOR SAND GRAIN 
BEDS (1-5 0) 
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So V 

.-lsl'.i.L0:'.3-ip of p.^vam-tcr- ’,.ith rnu^^‘'.-!^s3 c..'.c iit- 

> 

r-'.tio’.; TvYid noTiJ.iif •'irT!it;v — f'~/<^ Ji j 

r- ^/*^50 dlicused 

in do'call in tbe follov.l.ig su bn a 3 1 ion's. 

of -I' ifora 

Afe^ j-buted 

( i) ■-•'ifeirads:?’ -j "quiv.ai r'it s;n1 fazi'z rout^b’^-ss (Kg): 

xh:s i/Otal aif-'ot of roughn-'-iS of •' surfaa ' 0 3.n 
0.- rcor-T'^-j.-'t-d by a si'^gli rougbn^'cin sc-J.-' Kg. fbo 
roughn-03 sc.:do pax-.mzU:: K^/d^Q I- clotted against /. in 
:?ig. 10, Sho r. suits I'^port'd by Sculichting (1936), 
O’Loughlin 3t.;d. (1964), Dr/id (1980) and S^rln (1980) 
hcv.- ;.'lao baon plutt'd. ?or all theso c sna, lu is 3< en 
tb.?t th-. value of I’g/d^Q i icraasnr. with /, , attains 
a mvjcinius v.-lui- at - ,25 and then deer 'as as with 
further increase in >. . ^evid and Sai’in h:.vp. shown th 
foie h <0.1, Eg/d^Q Varies lln;,civ with , the 
constant o'. p:r-oportic«c0.ity le function of d^Q , 

.for: valu-'s /% X).1, Toaghnneg sc decra^ap*^. !?h? 
functional forta of Eg/d^o with A m.x7 b'. wrlet^n as 



Ks/dso^d-'/A) 



FIG. 10 EFFECT OF ROUGHNESS CONCENTRATION ON Kg 



3 ^ 






3/2 


(15) 


• 5 lo ■ r c* 


I i 


If "" 7? ’ : ' 


•y ' * * '^■*'■''3 ri~h ';xp‘'rin-’’t ’. 

f irlT W-- -I . -? .Uiia tod i- -jO. 


( tl) 3hif; ii .glocity sc -l^ 

!I?b' 3roaghTip.f!8 sc.ai'* r‘-pr'=i^''V''d =",3 1-^t < 


T loci '-7 3C 1- 


Au 


i3 furctJloc of r, ' of 

T* . den 

rou~bn?s3 surface (x ) ad flow prop'-:.:!.^ 

V 

3 vld 3nr* 3 rin sboncd fcb jt for oor''!: at: au.-f c=> g on-trr (a)» 
incroas s rin'"'"TT? with in " 


H->ao' , 'ih-' funotion -lyi - 


is HU^^OO^«»d *50 

b- f.irotic-i of tb" g^o-r-try. fblf? fu.-otlo- is plo'-ad 

^iu 3 G X in yiy. 11 . It mi 7 actad th V — .-,1 In ^ll 

V 

dr-cr^'i39a on Pith^r sldo of X =• •S29, alnll^r to t^c 
ViTlrtio-. of Eg/d^ with X • position of oopanm-'c-' ©f 

iu 

V'-.la--r3 ooincid* for all otsna wish X » 0.»CH% how^y"", 

■ tb-r-' n ■goiiajf 0 - of p^elt differ in o^ph s**, I; 

.n'’.7.h-’. a>w?rT«^ tb.'.t na^ltodow of p- >Jc.T'ilaoe 

grv3|l*:9M 

.'"M A’- ' ^ '*» , 



1 D Sand grain beds(Sand at nonunitorm 
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FIG.1! EFFECT OF ROUGHNESS CONCENTRATJON ON LAW OF WALL INTERCEPT 
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acatter. Th- differptics in ppik valu?s of glass b-ads 

(David ond Scb.lioh'tiQtJ data) is dus to regularity in 

shape and arrangement pattern. The difference in peak 

values of 3D gliss heads dataSa'''in, is indicative of 

regularity in shap.e and rardomn^ss in arraege'n''nt pattern. 

The present results adds th* findings of 3arln ch-it 

Irregulcrity in shape and arrangement pattern decreases the 

roughness scale - 1-ln ^ 0 . 7 * 

V >’ 


( iii) Theore tical h ed level S/dc r> 

Theoretical htd leVel (S) is referred to as th 
location of the apparent origin for velocity distribution 
abova th'' smooth surf-^v of the flat bed. 

The plot of s/d^Q against )\ is shown In Fig. 12, 

It is observJ;d that thr theoretical bod I'^vcl shifts 
from the top of th'. grains twoerds its bottom as the 
roughness concentration of th' bed dicreascs for 3DS 
bc'd ssries. For 0,925 D3 and 1.5 DS sand bfd 33ries» 3/^50 • 

is greater than on?, ihdicating that tfc- thoorrtioal bed 

% 

level lies above d^o level. This is du' to predominanoa 
of viscosity effect in these sand h;d sorlas. 



= 150 to 300 



40 


( no nuniformity In grain alas on the Rougtingsp of 
dgnaa ly covrire d by tb grains 

(i) Fikuradsr.fs roughnoss funotlon Kg/d^Q 

rikuradsv '»3 roughness function Kg/d^^ la plotted 
against nonuniformity par^rartcr ^'/d^Q for data of 1.5 B 
along with those of P,S. Mittal, 3. Mittal and Bavid in 
Fig. 13. It is ohservsd that as nonnniformity 'T/dijQ 
increases, roughness fhnetion Eg/d^Q also incrcasea, attains 
a peak value and then dooroases with further Increase in 
The value of at which maxim ;bi 


occurs varies with th: rmyc of 


^*^50 


for 




K. 


s 

>^0, has a unique curve for '^'“/d^Q, wheroaw 


^50 d, 


it gives different curves for 


50 




< 200. The slope 


of th' rising limhs of thes' curves gradually decrens- s with 

dj^ 

^ and that of the receding llnibs 


the decro.ase in 




remans constant. The peak values of all the bed series 
nre fairly constant. The present exporlmeatal results 
covers the wide gcq> left htetwooa th' datas of Bavld. and 
P.K. Mittal and S. Mittal i 


The shift in carves of tg/%o <^/«50 *®*’ 

different d^pt may ho attniwited to obange In th« 
of flow fro® roui^ turbulont non to Mtooth t«xWn«t tlw. 



3 D(S. MITTAL) 

4D,6D,8D {RK. MITTAL) 



EFFgCT OF NONUNIFORMITY IN GRAIN SIZE ON EFFECTIVE ROUGHNES 
SCALE K« 
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Zov; 

V'r, 1" ir.T 

b not'd :h t 

20 } 

r'i? r: -’'o' 

> 200, 

:zp ■ 

rim r.t.n It 1 

- in tt T., /io 

Zor 

c < '^50 '» 

' ^ V- 

< 70 .n t 


ir'ido , li tv. 

pr* ^ rt c v.n 


< 200 . 

y 


^50% 

^ m m mmrn 

. :’. r: r 
r • 1 on 


< -0 


cm ’ l(nr*.i3 .'n 
. It: ’*:glon 
:c 0 il ’ X'' 


( ii) scale 

T* 


ih’’ ^ foot nf roug^.7t=J3S c?»n I'^o b* •?XDJ*9S3'>d tV- 
AU 


vertic! 1 3Viift. ( ) in t’" v->loclt 7 fiintrLbatlon froraaoKjotli 


wa.i 30 V ;:na h r.all. velocity s''l.f': is plotted 

• 4 . %0 


^ .Tor noHunl-fom send ^?r~ln beds using 
pr s=‘'it experim ntal r->siilts ?long ’fti .'h tb** data of TJartl 
■ 3. 'littol IS shown in Fig. 14. It aiey be obnsrr -d 
tbnt data of certaia rnaee of nouu ilfomjlty follow 




a p'-jTtioular ouj-ve. !!!he r''glon ' ebc'rc >300, end 

V 

region 20 < d^Q^- < 10 follow a straight line ▼ rlatloa 

on semilog plot with slope equal to t/'«, In region 
d^Q+ <20, curves daap down, Infloatirg towards the region 
of s;i :)th turbul'^nt .flow. In th'^ r*^glon lD<-~ * <300 
tho curv'-'S crosses a.aoh other «»ept for eunr* baring low 
valu-a o-’’ noaaniforwt.tgr ( ' 9a^ preireat 

exp'^rlm3nb4i' inresttgritlon-'^i'e*® oi(^ out sprol’tllF to 
verl-?y wh-^ther such ..cKossd^ exists or not. 1&qperl?«'*at'a 






ria. 14 VARIATION OF A«/V* WITH FOR SetECTEO VALUES OF 

NONUNlFCmMITY PARAMETER <r^/4sQ FOR ALL BED SERIES 
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r^3ul':s d? prov= 
’Uhls ph®noa-na i 
due to visooslt/ 


exieSou3e of crosav.ig of v» curves, 
o -uc ■’o shift of rou^'H^ss cunrss 

. Hogubnr sc l-* -models wmtfi «i*r'*lop'd 


s utis ecu ' n tly . 


Thp function el fo-.-n - I In 

*50 ’■ ^ 


Is nupposod 


’GO be independent of —— — * and function® of ^ /d^ 
only. In order to study this Taii*tion 1 In d^. 

is plo'jtac' against 6’/d^ as s' own In the ?lg. 15, It is 


observed, tint .as nonuniforciity 6'/d^Q Incr’aaea, 

1 , - + 

increaseg, reaches a pe k ▼ due and 


Au"*" - d In dj-n 
oc 50 


then decreases with furcher inorease In cf/d^, ::jh«. 

value of <5yd^Q at which nnxi;au8i Au^ *• •^In d^ occurs 

varies wi^:h tbe range of 2i , for — J-S > aOO, 

, ' V V * 

"t* ^ 

"-x '^50 ^ uniba® curve fox 'S/d^, where ao 


it given di 'ferent curves for 


%0 


K 050 # yhp' p^%k 


+ 1 4 * 

of *--111 d^Q mm^lm smf for '^seb 

^50 

range o:. « 


The present '-rxpertii'eiit^” result cover®' the wide gap 
between the d-at; of F*X, ^ttal, Si^'^’flttal tiiid IKvld* 

Exist aac® of ^^ml^c'-ipr## .for d^ _ >200 war t>« 
rtttrlbut'^d'tto.witi#.''tiip^ictt JS-'tlf. flim «h»ar« ftpo« 





w% 







0-92 5 D 
0-39 D 
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TO.fS EFFECT ^ NONUNtFORMITY IN GRAIN SIZE OH THE LAW OF THE WALL 
WTERCEFT 
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rou ’n tur’Dulent to smooth tur'nule .i; flow, th-= curves 

snift wiGh p?jrti.oul&r values of Grain Sfc*'<x ReTaolds 
number. 

( iii) beg l -:vtl 3 /d g^ 

Theoretical beg level *£» I 3 mersur* 'i from the 
flat aurf.iCG on which the sand grains are stuck. The 
magnitude, S/d^^ is found to very in particular way for 

^*^50 

given -rp — » as shown in Pig. 16, The value of 

is maximum, for 6 /d^ = 0,05 and decreases with 

increase in C/d^. Ssperiments reported in literature 

also indicate the vilua of 2/d^ - 0.75. However, in the 

present investigation and invest ig'-»t ion carried out by 

David, 3 was found to lie above dcQ In aany cases. 

This eff'^ct may be attributed low values of * , 

' V 

or transition from rough turbulent flow to snooty turbulent 
flow, 

\ 

D, M odel f or Roughness Heals 

Prom the study mode In prevlons p«ragra|^» It was 
shown that roughness scales represented In the forn of 
^s +1 *50^* 

3 — or An - In Is fnactltm of pasraostsrs 

^50 X V 

representing 0»m0trf of %W bed lite wwtllmbsi 
concent rat ion, > Slid nonnnlforaiity psaroneter ^ ini 

' " ' • ■ ' . 4 . I 

■ ■ ■■ ■■■■ ■ ■ ' ^ 




F^t6 EFFECT OF NONUNIFORMITY IN ORAIN SIZE ON THE THEORETICAL 
8EO LEVEL 
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r'preo: r,tlng th. 
Raynolds nud^'r, , 


stat;; rioi^ n^d lv (Sr *c Sh*i<9Lr 
imod-il for i'oug’in 3 p.c-lf ‘9 


in t^r-rds of abovr p Tarn • tors h-s b'o.i at- v - 1-' •'* U tW 


following p.'Tt^Taphs. 


^ Modr;l bas'o d on th' g-pta otry o f the x^ouch bad profile 

Rrou' the literature study it was notr.d that a 

porticulax raprrsentative size for beds hating uoauGlforti 

s^d grain densely packed was shoosen* !Po nam- th‘ 

few such sizes ar-^ d^, dg^ and d^Q. It was shown by 

David, P.!K. Mitt cl And S. Mittal that along wltti 

poTtlculftr size, it is necessary to use seoond statistical 

peronetcr aa^i''ly standard deviation. Standard davl-tion *6* ' 
was 

/found^ to give betti-r represent at Ivb gr'oaciaaiis aniilii* 

Standard devi 'tion donot eonalder the perslstanC' 
occurrenoe of particular slz’. In order to inoliido this 
ch-ract eristic, the antooorrclatlon eai? powor apr.ctral stady 
of th send gi'aliB protrusions is Btade# IQrploal plots 
of autocorrelation and pow.-!r spectral censlti’ shown 
in Mg. 17. • It msy-bo obdervad tern th« |iito«iormlatloii 
curve that the slope of th > carve hefom sero orossing Is 
function of ^osBctry oT, send bed profile* ?tor ualfteH gr'ilnOf 
densely packed where ^rotruclons ir® aw ill aad olosoXy 
spaced, the sl'^S’ o^^'-^tho ' ^firtooojsrol^ition owrvs 



Autocorrelation function 



d ia 16 

OistofiCe XCmifil 


20 I 


3 5 7 ^ 


FIG. 17 AUTOCORREtATION ^ GRAIN HEIGHT 

■ . FOR ^5 P5 ■.■AliG■:^5 - 
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stp.^p. ^ ‘ protrusioris liz? and sp'^clng incr^as", 

-h slo'i o_ aai:ocorr‘^l . tlon curr" decr^as^'O# tbt 

•slop? of autocori-il.tion ca- - I 3 conalder-d to '•^^pr^'^ont 
bb3 spacii?: choraot'irigtic of s/n^ ilaa o-i sand b»d. 
la ord^T to utilise this iaforaotloa th r"pr«*«?'»ni:ition of 
esomearic so.il-=, th? slop’^ of aatocorr^l .^ion curve d»aots4 

as 'Sa» W is multipli'^d with th’ standard d--vlatlon »'r» and 

3 

used in fch? study. Table 2 shows the v-duea of 

for dlffv=5r’nc <S'/d^Q. From th^ study of th-- table 

5'h^a 

it m-w noted tb.;t ^eemnlns fnlrly const^mt 

for all th’ ^ gl'^n bed series* 'i?bls 

oonsta»oy character doiiot jrepresoit the eh uractsrlstlea 
variation of rougnaaas scalas* Hone? this apfuroacb Is not 
considered for further analysis* 


Thr sand surface is assutaed to i'eprP8‘'nt oo*bln«tlon 
of number of irregal r waves of eertaln h»lj^ end spaeing 
Th'- significant wwe hot^t, em be oonsldartd to represent 
tb' one of th’ g-'cractric^ prop«rtl«*8 of l.rspil--r wares* 
This signific^mt ware hei^ Hg' , «ay b? repr*-S'.it*jd 
according to Gartwrt^t ( 1 % 5 ) ss 
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TA 'LS 2 : iIJ!POCORHSL45fIOK SLOPl YAI«ES 


31) 

.925 D 

cr~ 

^50 

Sa 


cr* 

*50 


s. 

^50 

.051 

5.0 

1.50 

.047 

2.5 

1.04 

.0804 

4.25 

1.40 

.203 

1.92 

0.805 

.0899 

4.15 

1.40 

*642 

1.47 

0.85 

.191 

4.0 

1.60 

1.302 

1.92 

0.91 

.339 ' 

3.22 

1.40 

2.37 

0,00 

0.92 

.781 

2,62 

1.13 

5.05 

1.22 

1.51 

1 .317 

1.82 

1.00 

12.0 

2.63 

2.55 

1 .786 

1.42 

1.30 


0.591) 


3.26 

0.87 

1.30 

.052 

3.57 

1.47 


1.55D 


.257 

1.40 

0.67 

.05 

3.24 

1.05 

1.227 

1.30 

1.74 

.10 

3.0 

i.ao 

3.593 

2.0 

3.26 

.202 

2.2 

.90 

5.629 

1.25 

2.19 

. .35 

0,635 

.95 

J.56t» 

■ . , 

1.17 

.65 

0.67 

1.12 


0.1431 

f 

1 .25 

0.475 


44' 

1.61 


• 


.052 

.71 

1.75 

0.396 







.aoB 

1.79 

.82 

4.80 

0.179 






.vt- * ¥, . v - . ' 



34*48 





2^940 

9*18 



"is ; * ‘y 
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In i?Vt<i.cno is c’-Endaird dBvi^'ioi ^ as 

a = ( 1 . il )V2 

wh^re Sq, S 2 and S^ -.re z'rotb, s“oond ind fouirth spectral 

moments ^^^hich in a genercd form can be repres«nt<*d is follows: 
a 

Sn = S(a>) dCu (17) 

o 

where S(t-J> ) is the significant spectral density corresponding 
to frequency^ (in radians per second). Rie ralnes of « 
computed by this formula «?• talralated In appendix. 

Using this computed values of 0 ^ and a, the slgnlf leant 
.height of sand surface protrusions was* c<M^i>ated and la plotted 
in Fig. 18. It may be observed that this slgnlflosiit tislght 
in the form Hg/d^Q increases with InorewKS In attain# 

a maximum value then decreases with fhrtber Increase 
in cT/d^o* variation is alrallar to the variation of 

Kg/d^o 'Wi'th /^o* tariatlim of »»* 

cr/a^Q in the Flg^ 18 . Itonunlformltv/paireiiwt^ }mam m 
protruatons height namely Increases 'With 

nonuniformity Cf* <K>a^iln v*J.tte 4 bS 

then aft-'T starts dsai^-^iug with, licsTJa**** lnor/4g;0'* fIlIJi 
variation Is ~*o% o'b^.e5^.;(d. ;th^ of bads- fflio 



SYMBOL BED SERIES 


-O 


-O 


6 



o -o 


4> 


o 



“O 


Q Q Q Q 
LO 

cs« cn sr 
If) cr^ CO «- 

« ^ M * 

^ o o o 


6-0 9 


-^5 


#*o 
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mean height of protrutiQiB Increase with Inciiease on 

certain maximum xalue and then deomases 
with further Increase in This Tarlatlon U consla- 

t.uitiy observed for all s\nd bed aeries. 

2h*om these geometric -J. model study, It may he 
observed thrt the state of flow has to he t-^en into consider** 
at ion. in analysis. With this s» aim, the following a»dcls 
are davlopad for uniform grain randomly spaonA and for 
nonuniform gr.;ins densely pacted. 

(h) Model for Houghness of nniiora grains randomly nyaoad 

Roughness of uniform sand grains rMidonlT sp.used was . 
'-experimentally inv^stigat-d hy tavld for d^ <* 5*0 nm. Sarin for 
a.^Q = 0.925 mm and d^Q = 1.5 m® during present Inwestlgation. 

In all these inv- stigations the fallowing oonwm fm&tnres arw 
ohserv’^d. The maximum value of roughness soale lilce 
occurs \t or unround roughness conoeBtrxtloto 7i » 0.26 »aid slth'f'r 
side this value of , » _■ '^^^di^j'.deciPaases* 8 i 3 ®sitnde of 
the maximum value deortia^ea'-.wiltl* decm-3ai^ in ffitnln SSb&dSS 
Rejmoids number. . ifedelllng rirnglussn Wile 

is proposed in the following ^aph whloh inoorporatoe 
above features. »D^#M|ig if by xelKtIng theoretlosa 

bed • S'/d^o with. , */%o ’■itli 

^ and . 
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or.;tieril bed Ia , 

/ 50 'igilasrt it>ughnns8 

conc'^Ttration 7 s as shown In Fie, 12 «>« «>«. 

jjit,. K, I- wjs ol^erT'-d th^t 

■ih-. .ocpcnmcni; .! vlu^s of 

p-.^-.llol ourT>5 for ;-oh r,-ai-0 0I3, 3.0 m, 1.5 «, 

'^-nd 0.925 nim. Sbos^.- curv’a .... « 

Y a ^ Gr'iln 

Shear Reynolds nunbsr * 
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3mplrio.al equation given below 


^ ^ shown in Fig, 12 using an 


* 

^0 


386 


( Ji.. )“-♦ 


(f5^)1.13' m 2 

'7 


( 18 > 


I’he roughn ss scale: ropi-esentea In te-rms of shift in 

velocity ind theoretic l b-'d lavcl oro found to be r-el -tod 
dsrn F* 

wihh > -nd — os 


£jU 


"i- 


+ ■ 

wh^'r- 3 = :ind : 


In 8''".= 1,5 Bxp. 0,00535(a0f5 )^^^Cl10*-%) 

in) 


-V 


Here the p-oraTiGtor n*^oii^.iits tbe^ son of wmm hr^lf^t 
'".nd stand-.ird deviation of hoists as ' • 


t = ( 1 -^^) ( 1 .^) 1 :^ 


(r£ ^enV, 


50 ^50 


( 20 ) 


Using thn above fo^tioiis# i*© xnlatton fof sldfl l» tlia 

Tclocity scale is wnittoB.'-^'-v ' ■ . 


■ ' I' ’■ .t '/.U** -'V ‘ ‘ j" •' • . '. "js ‘ '*14^ t i* ’ 
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TAS.LE 3 : OOMFJTBl R0TJQISES3 33411 ?4HJBS 



5 DS 


1. 

5 D8 

n 0.925 US 

'"1 

Kg/d 

50 

xl 

Kg 


A 

V«so 

Sxptt, 

CJoaputsd 

Escptt 

. G!or:puted 


Oominited 

.01 

.28 

.3744 

.05 

1.05 

1.0973 

.10 

1.1 

1.66 

.02 

.55 

.6588 

.15 

2.0 

2-29 

.20 

1.85 

2.32 

.065 

1.70 

1-4735 

.26 

2.40 

2.84 

.30 

2.05 

2#66 

.150 

2.45 

2.2500 

.40 

2.33 

2.77 

.39 

1.00 

2.645 

.262 

2.95 

2.5100 

.70 

1.72 

1.47 

.47 

1.65 

2.223 

.394 

2.70 

2.9680 


■ 

0, 

.75 

1.25 

2.210 


2.2400 


.608 2.30 
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= ^ In : 


' ; 386 

,1/3, 


+ 1.5- Bsp(0. 00555(20+^) ^^^(110..^^ )) - 


The roughness scale Sg/d^ la related with aJbore 
equation using tte equations 9 and 21 as 




3^ ^ 3.215 I 366 1 Jl.5-*i5,.00535{a>*-'q) 

So (45^)^-’’ 


Tjaing the average d^o^ values for these three hed eerlee 
the val’ue of Eg/d^Q is calculated and teliol«iilie3U>iig with 
experimen+al values in table 5 shown helow. Ifhe 0(wa|>Uited 
values of Kg/d^Q are consistently higher la c«R5»arl»oB to 
experimental values. The niaxlwuffl value of oocare 

around ^ =0.3. 



c . Model for Roughness of sand beds having ooimalfO^ eaad 
grains densely papteed 


Experimental results of series of sand hods with U5 
as median diameter of the pres®*^ woritii. «H^^,,tlie -di^e Of. .. 
P.K. Mittal, S. Mittal and David for sand bed series for 
different laediaa grain di^etears reiigiiig S*P 9m 

'to 0.14. iw are' iise4 .in..t|>e of .jBOdel fOr e<aid 




1 
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beds with ncmmiform sand grains densely pacl^* froa the 
expGri?a'."!t3l data, referrln<! to Tig. 13, the following 
observations n?v7 be noted. The s.agnltude of 
increases rith inersase in ^/d^O, attains a aaxlaun walue 
and with fu-^ber increase m It decreases. TbM 

position at which Tiaxiiaunj value occurs appears to be sasie 


for higher values of grain shear Reynolds nuiiber 
C — > 200) and shifts with decrease in 




for 


m 

< 200. The magnitude of maxiiBum value of Z^/d^O 


remains fairlv constant in all the cases* With these as 
noteablG features, the model has been developed on stnllacr 
lines for uniform sand grains raadosly spaosd osse« 


•( i ) functional relationship hetweon theoretioal bed level 


dcQV, 


and ^50^* and 


Theoretical bed level g, mummtA from tt» flat! 
surface on which sand grains are stnoR is Amotion of 
geometry of thr saaad surface nMcly nod state of 

flow represented in terms of • |!roa tho experimental 

obsejnratlons It may he noted that aet portloul^ir valme of 
^50^* , the S/d^Q is naxiipim wh^ or* is sdnimwiy la 

pi^sent investigatt^ ■issmnimi nl.'lli 

increase in nommlfbiaiity Is olnMBrvoft 

in all cases. Hth en. nsssisp'tlixn' that' Is jntiNilj 
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f-3;:cT-an of crVd^Q for bighci* v-ilues cf ^50^* and 
f -nTblorr of bot- Mb fir So''« ^^00 

an '-.‘ripirioal ocuatlon Is propo'^od as shoim In Pig; 19 


.1/3 


•50 


= a-rpy: ^0.18 +-^ J'^d50 - (l+^^/d^o) 

*v „ 

l 3 cp 0 . 0144 ( 47-^50 


(25) 


Gonstdering 3/d50 is function of and 

the plot of 3/d50 with + d50 was made as shown in 


the ]?is. 20. 


It may be seen that 8/450 decreases nonotonicallp 
with increase in q~- + 450 and their fhnotional relationship 
may be represent d in the tom of 



4.4 

1 + 0 . 026 (^ 5 ^ + d^ ) 


®2 


(a*) 


(ii) Functional relationship anosg m* $ t'*', d^^ moA 

Using the present erperloental rewiXts, along tiith 
results of pavid, P.I. Mittal and S. mttail, thl 6 plot of 

- L In e'^ with lfTd^ is aiado in fig. 21. Ow 
p araoterer ^ «as ohosen with an idea 'thad d^n 

average width of protrusi<«» aaad 

of protrusions, ^e prodact of I50 repmsoats ths ^ 
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VACATION OF Ay*-| In t* with s/ <T*di 
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effective yxe a. of t h“ protrusions faclnf^ the flow. Th« 
parameter represfnts with Grain 3hc«p Rnjnold* 

nurfosr 'D-iicd on the s qu.ire root of avt'ragf* effective area. 
T he fuin tlon.al relationship between Ixn* - -Jl In and 
nia^r va'itten as 




0.03 (35-l(5?^d5J )J 


( 25 ) 


Another functional relationship has been defvelopsd 
using ( cj^ + d^o ) as parasaeter. This functional fo»i Is 
represented as 




J In S'*" = 3 - 9,6 mxp. - ** ^2 


( 26 ) 


where y = ((7^+ d^ )(t - Ijq?, ( 


m 


CT" 


)) 


Both ab^'ve equations are used in the #rveli^«at of 
model for roughness sosdLea, 


(iii) Model for rouglmegs 8Cii3-® 


i ■ * ' 

Rou^iness scales 

^’‘***' **^^ ^ 850 


wm 


related with CT/djo ^d d^ the 1(|8« 23 or 8|« 19 or S6, 

Using equation an eafprwsion iPCtt* %0 ^ 

written as 
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A.u^ <3150 “ ^ 2 ! 


(27) 


Im '3>:pre3Sioa for JTikuradse's xoug^nftsfi acnl« VSo 
;ia 7 be vjritten using Sqns. 9 rjid 27 as 




+ 3.215 r^i 03 ? Igl lip.OCfi or Xf^ 

(SB) 


Jtaiong the combinations of &| or B 2 with cep fg, 

from the computed Values it was found that a oonhlnatlofi of 

Bo with B-i gives a better representative values In 
^ with 

comparision ^ experimental values, Jn expression for 
A u"^ - / In dcQ and Kg/a50 ooBblnatlon 

me. .7 be vTitten as 


A ^ In dc-n “ InF— — — , l '\ I 

^ X ^ ll+, 026 ((f^> « 5 $ ) J 






^n(2.718+(.002B .a^ |l-»- 0 


and 

K 


. 05 { 55 -^<r*«^)i 

( 39 ) 


/ 


s _ 
d 


3.3 


3.215 I *.4 fcpfia 2.718K><»a H* 

5 Q «50 1 +. 026 (<r +^' ) 

[i-Bqp . 05 C 5 I 5'4 ^ )*] 

liable 4 below shows t;he ooagnrfeeil *b1u^ of V «56 1«i 


corre^onding ©3 
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2 43LS 4 : OOT-tPITTin) YOTSS 0? RmORmiSS SCAI.1S 


^■•g/ a^Q 


0.925 » 
^%0 


.051 

1.75 

3.25 

.05 

1.65 

2.20 

.047 

1.55 

1.95 

.063 

2.05 

3.73 

.11 

2.45 

» 

3.20 

.203 

2.70 

3.70 

.08 

3.45 

4.30 

.202 

3.40 

4.56 

•642 

3.90 

5.65 

.09 

3.69 

4.58 

.35 

4.20 

5.60 

1.302 

4.60 

6.03 

.191 

4.90 

6.13 

.65 

4.65 

6.25 

2.37 

5.00 

5.33 

.339 

4.50 

6.50 

1.25 

4.25 

5.78 

5.049 

4t0 

3.59 

.781 

5.40 

5.53 

1.75 

3.70 

5.13 

12^ 



1.317 

2.90 

4 i 41 

4.80 

2.05 

2.74 

( j . 



1.786 

2,60 

3.72 

0.39 B 





3.26 

1.80 

2.48 

.052 

1.30 

1.87 




0.14 D 


.257' 

2.^ 

2.78 




•052 

1.10 

1.68 

1.227 . 

3.30 

4.29 




• 208 

1.50 

1.94 

3.593 ’ 

4-65 

4,e9 


: ’ ■ 


•766 

2.00 

2.34 

5.629 

4-15 

'4*36 




2.iao 

3.05 

2.67 

9.567 

. 3.35 .>3^5 

k ... Ji . 



, 




, 1 . 

-T 





r .. -K ~ 
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?rom table it nay be observed that position of peak 
valtins coincide with position of peak occurrence of the 
experimental results. The trend In the varlattcn follows the 
exp'-rimental values. However, the nagnltudes of peak values 
are higher in comparision to corresponding experin-ntal values, 
Also the magnitude of peak values increases with decrease in 
^50^* , reaches a highest magnitude around d^ * ** 100 am^ 
d^reases with further decrease in d^ . values- 


S. General Gomments 


Development of models based on gemetsric parasieteQrs 
like OE" Hs found to be inadequate to iw^resent 

the roughness scales. Along with these gsofflwrtrlo psiranstiTf 
the median size also has to be considered, 3Wrth®Pf oonai*" 
deration of the effect of state of flow baa to bo BaAo* 


Jrom the models developed for untfox« sand gi^n 
randomly spaced, and nonuniform sand grains danaalsr 
packed, the trend in variationsooiscias with «q»oriisafital 
values very nicely. The order Of magnltudS of Tm 0 m»m 
scales are also within t » per oont of tke Wr*iM»tai 
values. Howev r, a farther reli^iowt W 

help to achieve a close agres^o leodalb* 


An effort was to .^a. anwi^ fo*f 

UnifoCT r^dosily sapced- aad' fwdkad 
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based on the otrtlstical paraneters of aurfic-^ protrusions 
na^-ely, oi protrusion heights ana standard dOTlatton 

of protrij.sion hoighT:3. This was not very sncoossfUX and 
h^^noc not reported. 



chapter IV 


GOITCLl'SIOITS A!rD HSCC.'T'EIDjfflICN'S 

A. Q-en eral 

Roughness of sand beds was investigated by aeasuring 
velocity distribution on flat sand beds in wind tunnel. 

Sand beds used are of two types nanely, uniform sand grain 
size of 1.5 mm randomly spaced with 5 different rougbaess 
concentrations and beds having densely paclced nonuniforB 
sand grains wii-h median value of 1,5 hhb for 8 different 
nonuniformity values. Experiments were specially aimed to 
investigate the transitional state of flow. Flow parMetmps 
like shift in tlieoretic-il bod level E, shift in velocity 
scale u bed shear velocity of the flow, were cwapited 
from me.an velocity data in tho wall region using statlstloal 
method. The velocity distributions were malysed using 
these parameters and showed how rou^iness concentrations 
and nonuniformity parameters affect the velocity dlstrlbutjUKsyi. 
Velocity profiles are found to be represented by aoatles 
representing the roughness; in terms of lllDara&f#*s 
equivalent sand grain roughness Ig/d^Q and shift in velocity 
scale £h.u^ *-J~ In d^Q and shift in length scale S/dj^* 

These scales are functionally to O'/d^ or,^ md y. * t 
using the present experimental data along with the data of 
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I)avic>, .Tud d. -'Ittal. Iho follo'.vine ar- the 

ma,-c'r coac ions b^3yd r.n tV:e an il''-s ''.‘3 of tl*' e:-:pcrlment.al 
dab a • 


3 . 'donolosions 

1 . Ttic axistance of transition in velocity scale Au"^ 

in the region 70 < < 200 with different ^ /d^Q 

is established, experimentally ^ 

2. The velocity scales and length scales represented as 
“f" / 4“ 

Au - ^ In S and O/dcQ arc found to toe ftonct- 
ionally related vdth or d^dcQ and — . 

^ V 

5. The models for prodicting the roughness of sand 
beds in terms of shift in velocity scale 
and :Tikuradse’s equivalent sand grain roughaesa 
■a?.’e developed. These models predict the experinental 
results fairly well. 


0. Recommendations 

IJodels developed for roughness scales are sultatol® 
for two separate cases namely uniform sand grains raa^orty 

n 

sps>,c6d. ^.nd grains dsiissly A 

"bo be used for th#' both cases is needed to bs developed* 
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3.T.sed on tho some analrsis c-jrried out, it may be 
suggested thet models should be b-ised on the stgrtlatioal 
parameters representing the surface of the beds. These 
statistical parameters mar be of mean of heights and 
stcjidsxd deviations of heights of protrusions of roughness 

elements , 
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Bed Series 

Ib 

^50 

Ha 

r* 

•^0 

SgxlO^ 

S^xlO® 

m 


1 

2 

3 

4 

5 

6 

7 

3.0 BS 

0.01 

0.56 

0.50 

0.30 

2.05 

0.25 

0.67 


0.02 

0.77 

0.71 

0.40 

1.41 

0.09 

0.67 


0.€?«5 

1*.08 

1.57 

0.96 

4-49 

0.04 

0.69 


0.15 

1,29 

1.88 

1,07 

2.24 

0.10 

0.71 


0.202 

1.57 

1.88 

1.58 

2.70 

0,12 

0.79 


0.594 

1.62 

1.67 

1.74 

2.81 

0.13 

0.80 


0.608 

1.56 

1.47 

2.02 

0.86 

0.02 

0.90 

1.50 BS 

0.050 

1.20 

0.88 

0.25 

5.32 

24.11 

0.99 


0.150 

1.05 

1.90 

0.21 

0.73 

0.05 

0.72 


0.260 

1.25 

1.65 

0.25 

0. 35 

0.01 

0.75 


0.400 

1.57 

1.46 

0,38 

0,46 

0.01 

0.79 


0.700 

1.66 

1,04 

0.31 

0.C9 

0.001 

0^86 

0.925.BS 

0.100 

0.85 

1.29 

10.10 

83.70 

12.35 

0.65 


0.200 

t.01 

1.35 

18.16 

148.01 

20.99 

0.64 


0.300 

1.12 

1.83 

8.67 

22.77 

1.06 

0.66 


0.390 

1.14 

1,58 

15.^ 

5r.62 

1.54 

0,65 


0.470 

1.14 

1.48 

8.54 

7,41 

0.10 

0.60 


C »150 

1.23 

1.02 

12.17 

10.99 

0.17 

0.64 
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Bed Series 

1 

2 

3 

4 

5 

6 c 

x10> 

7 

5-00 D 

0.051 

1.40 

1.53 

1.21 

o.oe 

0.51 

0.95 


0.080 

1.03 

1.99 

4.54 ’ 

' 1.05 

6.34 

0.79 


0.090 

1.04 

3.31 

4-54 

1,17 

7.44 

0.77 


0.191 

1 . 26 

2.86 

4.59 

1.03 

6.10 

0.71 


0.559 

1-5? 

1.92 

4.71 

1.07 

6.21 

0.78 


0.781 

1.37 

2.48 

4.67 

1.04 

6.21 

0.79 


1-317 

1-73 

1.68 

4.90 

1.05 

6.04 

0.79 


1.786 

2.96 

0.87 

5.44 

1.07 

6.15 

0.81 


3.360 

5.43 

0.33 

12.68 

1.24 

6.65 

0.^ 

1.5 3- 

0.050 

1.26 

1.31 

0.50 

0.012 

0.C75 

0.98 


0.110 

1.13 

2. 17 

3.46 

1.160 

7.f10 

n.7l 


0.302 

1.17 

2.91 

3.49 

1.150 

7.69 

0.715 


0.350 

4.27 

0.98' 

5.01 

■,1.58 

10.50 

0.725 


0.,650 

, 4.83, 

0.96 

5.36 

1.65 

1#.60 ■ 

0.723 


1.250 

^ 6.90 

0.61 

6. 14 

1,85 

12.35 

0.750 


1.750 

8.50 

0.50 

13.60 

1.95 

11,90 

0.790 


4.SOO 

16.95 

0.12 

21,48 

2.32 

12.82 

0.890 
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Bed Series 1 

2 

3 

4 

5 

6 

7 

0.59 B ' 0.052 

1.53 

0.85 

1,66 

0.043 

0.09 

0.933 

0-257 

1.89 

1.16 

3.00 

0.090 

0.055 

0.974 

1.227 

5.28 

0.625 

18.87 

2.110 

103.890 

0.990 

3.595 

5.51 

0.84 

16,61 

2.682 

14.960 

0.843 

5.629 

5.93 

0.70 

15.58 

2.480 

14.500 

0.854 

9.567 

4.13 

0.80 

5.^71 

1,200 

6.760 

0.790 

0.14 B 0.052 

1.712 

0.64 

0.46 

0.010 

0.039 

0.970 

0.208 

1.330 

1.12 

3.37 

^.010 

6.080 

0.703 

0.766 

3.205 

0.624 

3.917 

1.170 

7.230 

0.718 

2.120 

3.700 

0.823 

4.219 

1.259 

7.889 

0.723 
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Bed Series 

1 

2 

3 

4 

5 

6 

7 

1.5 BS 

0.05 

0.163 

0.476 

1,531 

0.121 

19.55 

6.24 


0.15 

0.245 

0.566 

1.521 

0.146 

17.35 

7.12 


0.26 

0.416 

0.687 

1.520 

0.149 

10.38 

7.65 


0.40 

0.687 

0.770 

1.530 

0.124 

7-42 

4*05 


0.70 

0.960 

0.752 

1.518 

0,155 

.6.92 

3.83 

0.925 BS 

0.10 

O.IH) 

0. 326 

0.986 

0.082 

17.06 

11.54 


0.20 

0.194 

0.395 

0.992 

0.065 

10.47 

4.66 


0,50 

0.234 

0.4?4 

0.967 

0.080 

8.84 

4- 44 


0,39 

0.2^ 

0.450 

0.969 

0.070 


1.91 


0.47 

0.3U 

0.464 

0.964 

0. 0 B 8 

6.26 

1.32 


0.75 

0.352 

0.480 

0.986 

0.083 

5,77 

0.636 


Bed 

Serif's 


2 

3 

4 

5 

6 

7 

1.5 

.9 

0.05 

1 .66 

0.505 

0.794 

0.554 

7.18 

5.70 



0.11 

1.68 

0.618 

m 04^ 

0.417 

6.71 

2.56 



0.20 

1.75 

0.630 

0.992 

0.715 

6.89 

2.88 



0.35 

2.12 

2.510 

1.940 

3.195 

5.75 

2.28 



0.65 

2.35 

2.590 

4.172 

8.950 

7.75 

3.23 



1.25 

2.93 

3.570 

4.585 

4.929 

7.55 

5.60 



1-75 

5-25 

4.250 

6.750 

6.870 

7.03 

3.12 



4.» 

8.50 

7-324 

9.562 

8.170 

6.67 

2.78 


